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ABSTRACT 


The evolution processes of double layers have been studied in a series of laboratory experiments. It was 
found that the existence of virtual cathode-type potential wells at the electron injection boundary was the dominant 
triggering mechanism. The rapid growth of the potential well led to collisionless ion trapping and the establishment 
of the necessary trapped ion population. For double layers with small potential drops, collisionless ion trapping 
actually induced ion-ion streaming instabilities and the formation of ion phase-space vortices. In this regime, the 
system often exhibited relaxation-type oscillations which corresponded to the disruption and the recovery of the 
double layers. 


I 


I. INTRODUCTION 


Much of our recent understanding of double layers has come from laboratory experiments and numerical 
simulations which had rather limited system dimensions. The system boundaries are often in close proximity with 
the double layer electric field, thus affecting almost all aspects of double layer physics. The situation is obviously 
different in space plasmas where boundaries are not well defined and often far away from the regions of possible 
double layer formation. In order to extrapolate the results from laboratory and computer experiments to the space 
context, it is important to understand the role of the system boundaries on the formation of double layers. 

Most double layers experiments (Quon and Wong, 1976; Leung et al., 1980; Singh and Schunk, 1983; 
lizuka et al., 1979; Saeki, et al., 1980) have utilized the injection of a drifting electron species to trigger the forma- 
tion process. It was found that a necessary condition for double layer formation is that the electron drift velocity v d 
exceed the thermal velocity v te of the ambient electrons. This condition results in the belief that the Buneman in- 
stability with an instability threshold of v d =’ v te was the triggering mechanism for double layers. However, double 
layers with potential drops <I> > T e /e, the electron temperature divided by the electron charge, have been observed 
(Hollenstein et al . , 1 980) experimentally with v d as small as 0.2 v te . Ion-acoustic turbulence instead of the Buneman 
instability was expected to be the triggering mechanism for double layer formation in that experiment. 

Numerical simulation (Sato and Okuda, 1980) of double layers with v d v te have found different results. No 
double layers with 4» > T e /e were found. Rather, a new class of double layers with non-montonic potential profiles 
and <}> T e /e was found. These double layers were always preceded by negative potential pulses and associated with 
current-driven, ion-acoustic turbulence. As such, these double layers have been identified as “ion acoustic” in order 
to distinguish them from the conventional double layers. Since an electron drift velocity of v d > v, e may not exist in 
space (e.g. , the auroral plasma), the ion-acoustic double layers have also become a subject of considerable interest. 

In this talk, we will discuss previous (Hershkowitzetal., 1981) as well as new experimental results in order 
to identify the formation mechanism of double layers in our triple plasma device. We begin with the roles of the 


PRECEDING PAGE BLANK NOT FILMED 


35 


boundaries on the steady state characteristics of double layers. It is shown that the drifting electrons provide the 
initial space charge for double layer formation, while the trapped ions determine the evolution process and the shape 
of the potential profile. It is also shown that the growth of virtual cathode-type potential wells at the electron injec- 
tion boundary is the triggering mechanism for double layers. Collisionless ion trapping by the potential well is 
found to be the main process for the establishment of the necessary trapped ion population. 

Furthermore, double layers with <J> T e /e are shown to be unstable to the evolution of ion phase-space 
vortices from ion-ion streaming instabilities. In this regime, the double layer potential profile has a strong resem- 
blance to ion-acoustic double layers. This result may represent an alternative explanation of the small electric field 
signatures observed in the auroral plasma. 


II. STEADY STATE EXPERIMENT 


It is useful to first describe the operational characteristics of the triple plasma device in order to get some 
insights into the sources of particles that support the double layer. The triple plasma device consists of two source 
plasmas bounding a target plasma. Each source is separated from the target chamber by two grids. Plasma potential 
in each chamber is determined by the bias voltages of the grid and the internal anode. The source plasmas arecreated 
by filament discharge in argon gas (operating pressure P c *£ 1 X 10 4 Torr) with density n s = 10 9 -*■ 1 0 10 cm 3 and T e 
~ 2 eV. The ionizing electrons are trapped by surface multidipole magnetic fields in the region closed to the fila- 
ments so that they cannot reach the target chamber and produce plasma there directly. We have also confined our 
study on double layers with <j> 10 V. These procedures ensure that ionization effects are minimized in our double 

layer experiment. A schematic of the triple plasma device is shown in Figure 1 . 

Stable double layers with cj> =S 5 T e /e can routinely be achieved using the boundary conditions shown in 
Figure 2. We chose to investigate these smaller double layers in order to limit the accelerated electron beam energy 
to below the ionization potential of argon. From the boundary conditions shown in Figure 2, we expect ions to only 
come from the high potential source. These ions are usually pre-accelerated into the target chamber by the potential 
difference between the high potential source plasma and the target plasma. These ions are further accelerated by the 
double layer into the low potential side. These beam ions either exit the target chamber at the left boundary or they 
charge exchange with neutrals and form cold ions (T, ~ 0.3 eV). The cold ions, once formed, are confined electro- 
statically by the potentials of grids B and C. Although the charge exchange reaction rate is relatively low at our 
operating neutral pressure, the cold ion density accumulates to a significant fraction of the beam ion density due to 
their long confinement times. 

Electrons which enter the target chamber from the high potential source are those in the tail of the Maxwel- 
lian source distribution function. These electrons have almost no drift energy, thus becoming the thermal electron 
species in the target plasma. This contrasts with the situation at the low potential end. There the tail of the source 
distribution function, which is energetic enough to get over the barrier provided by grid A, is accelerated into the 
target plasma. This results in an electron drift with the drift energy determined by the potential difference between 
grid A and the target plasma. 

The boundary conditions in this experiment are believed to play the following roles: 

1 . The high potential side boundary ensures that the ions will enter the target chamber with a flow velocity 
Uo > C s . This situation is quite similar to that of a sheath at a plasma boundary. Downward curvature of the plasma 
potential requires an ion drift velocity u e > C s . Since the high potential side electrons can be treated as approxi- 
mately isothermal, the “Bohm sheath criteria” applies in this case for the double layer as well. 
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2. Grid B acts as a potential barrier for the low potential source ions as well as for the charge exchange cold 

ions which formed between the double layer and grid B. Since the height of the potential barrier of grid B is roughly 
20 times the ion temperature (Tj = 0.3 eV), no low potential source ions are expected to enter the target chamber. On 
the other hand, the high potential source ions are at a much higher energy than grid B and can exit through grid B into 
the low potential source. As such the only source of thermal ions in the target plasma appears to be that of the charge 
exchange ions. & 

3. The low potential side boundary allows only electrons to drift in from the left. Such excess electron 
space charge may be neutralized only by the ion beam and the charge exchange ions. 

4. There is no externally applied electric field across the target plasma since grids B and C are at roughly the 

same potential. The formation of double layers is a result of the particle flow rather than that of an external electric 
field. 


Using the experimental boundary conditions and the particle distributions at the sources, it is possible to 
determine the potential profile across the target chamber by solving the Vlasov-Poisson equations. The details of 
such calculations have been described in an earlier paper (Hershkowitz et al., 1981) and will not be repeated here. 
Rather we will point out some results which are relevant to our present discussion. A typical solution of the target 
plasma potential profile and the boundary conditions employed is shown in Figure 3. The model has grid potentials 
similar to those shown in Figure 2. The double layer is formed in the region x L x Xr where <J>(x L ) = 0 and 4>(x R ) 
— <J>d- Using the dimensionless variable i|j = e<f>/T e , the density of the free ions n if and trapped electrons n e , entering 
from the high potential source are, respectively: 


n ir(W 4^-*) 


erfc 


L i 


( 1 ) 


and 


n et M = N e ( ^'^ 2) (1 - eifcy/f^) 


( 2 ) 


The density of the free electrons entering from the low potential source is: 


, , . _ N \p-\p . 
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where »Jj| , i|> 2 and v|i A are, respectively, 
lized to the electron temperature. Both 
N. 


the low potential source, high potential source, and grid A potentials norma- 
source particles are assumed to be Maxwellian distributed with equal density 
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The density of the charge exchange ions cannot be calculated from the boundary conditions, thus it can be 
treated as a variable or: 


n it OJ/) = 0n e (0) e 


-CTJT.# 


(4) 


and 


n. f (0) 

0 = 1 

n e (0) 


(5) 


P is a parameter which depends on the density ratio of the trapped ions to the beam ions and n e (0) is the total electron 
density at *]> = 0. 

We show the dependence of the double layer on the trapped ion density with 0 varying from 0 to 0.35 in 
Figure 4. As the trapped ion density increases, the double layer becomes more detached from the low potential side 
boundary. Since P = 0.35 corresponds closely to the potential profile in the experiment, it is possible that a signifi- 
cant amount of charge exchange ions are trapped by the double layer at the low potential side; i.e., p = 0.35 
corresponds to a trapped ion/beam ion density ratio of 54 percent. The trapped ions neutralize the excess negative 
space charge created by the drifting electrons, thus maintaining a uniform plasma potential at the low potential 
region of the double layer. 

The contributions of the various particle species on the double layer space charge are shown clearly in 
Figures 5a and 5b where the charge density profile and particle density profiles are plotted versus axial distance. As 
discussed earlier, the ion beam provides the positive charge density for the downward curvature of the double layer 
at the high potential side, while the drifting electrons supply the negative charge density for the upward curvature of 
the double layer at the low potential side. 


III. THE FORMATION MECHANISM 


In order to understand the triggering mechanism for the double layers in our experiment, we examine the 
temporal evolution of the target plasma potential profile with v D =£ v te and v D > v te . An extra grid is installed at the 
low potential side to facilitate the pulsing of the drifting electrons. The boundary conditions for this experiment are 
shown in Figure 6. A steady state target plasma with n e — 10 7 cm' 3 is extracted from the high side source, and the 
target plasma potential is quite uniform axially with 4 >t = 4 V. Low side source electrons and ions are normally 
excluded from the target plasma by the potential barriers of grid B (biased at -30 V) and grid C (biased at + 12 V), 
respectively. 
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At time t = 0, grid B is switched to the ground potential and the low side source electrons are accelerated 
into the target plasma by the potential difference between <|) T and ground, i.e., 



When the un-neutralized electron stream enters the target plasma, the entire target potential decreases rapidly from 4 
to 3 V in 10 |xs which results in v D = 1.2 v te . The temporal evolution of the target plasma potential profiles, as 
obtained with an emissive probe using Boxcar interferometer averaging technique, is shown in Figure 7. A potential 
well begins to form near the electron injection boundary at t = 50 pis. The potential well grows deeper and widens 
into a double layer at t > 400 pis. The amplitude of the double layer is <f> = 1.1 T e /e and appears to be quite stable. 
This result can be interpreted as follows. 

Electron injection from the low side source creates excess space charge at the injection boundary, and a 
virtual cathode-type potential well is formed to limit the injected current. The growth of the potential well is accom- 
panied by ion trapping in the potential well. As the density of the trapped ions increases, the double layer becomes 
detached from the electron injection boundary, in agreement with our earlier result on the effects of trapped ions (see 
Figure 3). Notice the double layer formation time of T DL 400 pis is considerably shorter than the charge exchange 
time of T cx = 1 ms in this experiment. At such, the trapped ion population cannot come entirely from the charge 
exchange ions which fall into the potential well. A more possible source is the neighborhood ions which fall into the 
well during its growing phase. These ions will actually get accelerated down the potential well with energies 
depending on their locations in the potential well. 

We further decrease the drift velocity of the injected electrons by decreasing the target plasma potential to cf> T 
— 3 V. When the un-neutralized electron stream enters the target plasma, <j> T decreases from 3 to 1 .5 V in 10 pis. As 
shown in Figure 8, a potential well is once again formed near the electron injection boundary at t = 50 pis when v D 
= 0.7 v te . At t > 75 pis, a small double layer with cj) — 0.5 T e /e has formed. However, in contrast with the earlier 
experiment, the double layer decays into an ion hole-like potential well. Note the similarity between the potential 
profile at t = 150 pis and an ion-acoustic double layer. 

The time history of the plasma potential (<J> L ) at an axial distance of x = 10 cm, the electron current flow 
across the target chamber from the low potential source (I c h) 5 and the ion saturation current (Ij L ) at x = 15 cm are 
shown in Figure 9 in order to illustrate the double layer formation processes. At t > 50 pis, cf) T = IV, and v D = 0.7 
v te , the growth of the potential well corresponds to the abrupt decrease of <j> L . On the other hand, I eH continues to 
increase due to the injected electron current until 4> L becomes negative where I eH begins to decrease rapidly. As <j> L 
reaches a minimum at -1.0 V, IeH returns almost to the level at t < 0. c|), subsequently becomes slightly more 
positive, and an intense low frequency noise appears in I iL which corresponds to the evolution of the ion hole-like 
pulse. 


A similar evolution process is observed when we increase v D just slightly. As shown in Figure 10, the double 
layer breaks into one or more ion hole-like pulses. The long time history of this experiment is shown in Figure 1 1 . 
The ion saturation current exhibits relaxation-type oscillations in time with a period roughly characterized by the 
transit time of the ion hole-like pulses across the target plasma. The relaxation oscillation corresponds to the evolu- 
tion of the double layer from virtual cathode potential well and the subsequent decay into ion hole-like pulses. When 
the pulses reach the high side boundary (e.g., the ion-hole velocity is the order of the ion thermal velocity), the 
process repeats itself. 
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The breaking of the double layer into the ion hole-like pulses can be understood as follows. The magnitude 
of the virtual cathode potential well <J> W is formed to limit the injected current. Since the potential well must become 
a potential barrier to the injected electrons in order to limit the current, 4> w = (vq/v^) 2 T e /e. The potential drop of the 
double layer 4> - <j> w ; also, we have <J> = (v D /v te ) 2 T e /e. 

As shown in Figure 9c, the injected current I eH and v D are reduced to very small values as a result of the 
formation of the double layer. We believe the growth of the potential well and the double layer formation also 
triggered bursts of counterstreaming ions which are accelerated down each side of the potential well with an average 
velocity: 



This results in a counterstreaming or “tuning fork” ion phase space configuration at the double layer front. 
As reported in many numerical and experimental studies (Pecseli and Trulsen, 1984; Chan et a!., 1984) of ion- 
acoustic shocks and ion holes, the ion-ion two-stream region becomes unstable when v b C s and evolves into one or 
more ion phase-space vortices. 

As we have observed double layer formation with v D =£ v Ie , it is doubtful that the Buneman instability plays 
any roles in triggering the formation of double layers in these experiments. When v D < v te , no steady double layer 
exists as a result of the ion two-stream instability and the evolution of ion phase-space vortices. 


IV. DISCUSSION 


We have reviewed results from a series of laboratory experiments concerning the formation of double layers 
in a triple plasma device. In steady state, the double layer electric field is sustained by the negative space charge of 
the drifting electrons and the positive space charge of the ion beam. The low potential boundary condition permits 
the injection of an un-neutralized electron species which space charge is crucial for the initiation of the virtual 
cathode potential well. The ion reflecting grid (grid B in Fig. 2) plays two roles; first, to prevent the low potential 
source ions from entering the double layer and second, to confine the charge exchange cold ions in the low potential 
side of the double layer. The charge exchange (trapped) ions are needed to neutralize part of the drifting electrons, 
thus allowing the double layer to move away from the low potential boundary. 

The formation phase of the double layers is associated with the growth of virtual cathode-type potential 
wells at the electron injection boundary. The formation of the virtual cathode potential well is a result of the lack of 
neutralizing ions at the electron injection boundary. As long as the injected electron density is sufficiently high, the 
potential well will form independent of v D and it need not be associated with instabilities. We have clearly shown the 
formation of double layers with v D < v te which is below the threshold of the Buneman instability. 

The movement of the double layer electric field away from the electron injection boundary is probably 
caused by the accumulation of the trapped ion density at the low potential side of the double layer as demonstrated 
by Figure 3. Since the double layer formation time is much shorter than the charge exchange time, the source of 
trapped ions is more likely coming from ions in the neighborhood of the potential well during the growth of the well. 
As such, these ions are accelerated down the potential well with a maximum velocity of: 
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(8) 



When the magnitude of the potential well <|> w < TVe, the ions become two stream unstable because v b < C . As 
such the double layer decays into ion phase space vortices with potential structures that resemble ion acoustic 
ouble layers. This situation is similar to the auroral plasma condition where small electric signatures (Temerin et 
al., 1982) are often observed along with counterstreaming ions. 

For the case of v D > v te , the depth of the potential well <j> w > T e /e and results in a stable double layer 
formation. For the stable double layers, the charge exchange ions will be the main fueling source for the trapped ion 
population in steady state. In that case, the main loss mechanism for the trapped ions is radial diffusion to the side 
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Figure 2. Axial potential profile of the laboratory double layer including the boundary condition. 


Figure 3. A solution of Poisson’s equation d 2 »J>/dx 2 = n ef {ij>) + n e ,(4») - n i |<ij»)-nj l (i|>) where *1/ 
ec{>/T e ; x = x/\ D ; n if , n et , n ef , and n u are given by equations (1) through (4), respectively. 
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Figure 8. Temporal evolution of the target plasma potential profile when v D > v te . 
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Figure 9. Time history of (a) the target plasma potential (<j> L ) as measured by an emissive probe at 
1 0 cm away from grid C; (b) the electron current (I eH ) collected by a single-sided Langmuir probe 
facing the low potential source; and (c) the ion saturation current (I iL ) collected by a wire probe at 1 5 
cm away from grid C. 
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Figure 10. Temporal evolution of the target plasma potential with v D < v te and the double layer 
decays into more than one ion hole-like pulses. 
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